Parametric Autoresonance in Faraday waves 
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We develop a theory of parametric excitation of weakly nonlinear standing gravity waves in a tank, 
which is under vertical vibrations with a slowly time-dependent ("chirped") vibration frequency. 
We show that, by using a negative chirp, one can excite a steadily growing wave via parametric 
autoresonance. The method of averaging is employed to derive the governing equations for the 
primary mode. These equations are solved analytically and numerically, for typical initial conditions, 
for both inviscid and weakly viscous fluids. It is shown that, when passing through resonance, 
capture into resonance always occurs when the chirp rate is sufficiently small. The critical chirp 
rate, above which breakdown of autoresonance occurs, is found for different initial conditions. The 
autoresonance excitation is expected to terminate at large amplitudes, when the underlying constant- 
frequency system ceases to exhibit a non-trivial stable fixed point. 

PACS numbers: 47.35.-K, 47.20.Ky, 05.45.-a 
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I. INTRODUCTION 

There are two ways of driving a classical nonlinear os- 
cillator by a small oscillating force: via either external or 
parametric resonance. In both cases, the initial growth of 
the amplitude of the oscillator is arrested, even without 
dissipation, when nonlinear effects come into play. This 
is due to the fact that the natural frequency of a non- 
linear oscillator is amplitude-dependent, so a mismatch 
between the (invariable) driving frequency and the nat- 
ural frequency appears 

To overcome the nonlinear mismatch and maintain 
phase locking between the driving force and the oscilla- 
tor, one can slowly vary the driving frequency with time 
so as to achieve a persistent growth of the oscillations. 
This simple and versatile method is called autoresonance: 
either external, or parametric. To emphasize the differ- 
ence between the two, let us watch a child on a swing. 
When a parent pushes the swing (once in each cycle) , he 
gradually increases the time interval between the pushes 
as the swing amplitude grows. Here he employs external 
autoresonance. On the contrary, when the child swings 
himself (he achieves it by moving the position of his cen- 
ter of mass up and down twice in each cycle), he gradually 
increases the period of these modulations as the swing 
amplitude grows. This is parametric autoresonance. 

The simple model of a nonlinear oscillator, excited via 
external autoresonance, has found numerous applications 
in physics, see e.g. Ref. 2j for a brief review. The ex- 
ternal autoresonance scheme has been also generalized to 
systems with an infinite number of degrees of freedom, 
such as nonlinear waves 0. LH Q and vortices §. In 
contrast to the external autoresonance, parametric au- 
toresonance has received much less attention 0. In this 
work we generalize to nonlinear waves the theory, devel- 
oped in Ref. for a nonlinear oscillator. Specifically, 
we show that the parametric autoresonance mechanism 
can be used for driving nonlinear standing gravity waves 
with a steadily growing amplitude on a free surface of a 
fluid. 

Michael Faraday || was the first to observe that, when 



a tank containing a fluid is periodically vibrated in the 
vertical direction, a standing wave pattern forms at the 
free surface of the fluid when the vibration frequency is 
twice the frequency of the surface vibrations. This phe- 
nomenon is a classic example of parametric resonance, 
because the vertical acceleration of the tank - an intrin- 
sic parameter of the system - depends on time via the 
periodic vibration. Lord Rayleigh |9j carried out a fur- 
ther series of experiments, which supported Faraday's ob- 
servations, and also developed a linear theory for these 
waves in terms of a linear Mathieu equation. Benjamin 
and Ursell 01 advanced the linear theory further. Sub- 
sequently, Miles mElIp, Douady [jj, Milner [Tp, 
and Decent and Craik [lfj formulated a weakly nonlin- 
ear theory based on amplitude expansion, while some of 
these and indeed numerous other works dealt with exper- 
imental studies of Faraday waves. 



Being interested in parametric autoresonance, we add 
a new dimension to the problem of Faraday waves and in- 
vestigate weakly nonlinear standing gravity waves formed 
when the vibration frequency is slowly decreased (chirped 
downwards) in time. We show that the negative fre- 
quency chirp causes a persistent growth of the wave am- 
plitude. Like in other instances of autoresonance, the ex- 
act form of the frequency chirp is unimportant once the 
chirp sign is correct, and the chirp rate is not too high. 
The autoresonance excitation is expected to terminate at 
large amplitudes, when an underlying dynamical system, 
corresponding to the case of a constant frequency, ceases 
to exhibit a non-trivial stable fixed point. 



Here is the layout of the rest of the paper. Section II 
presents a brief overview of theory of weakly nonlinear 
Faraday waves with a constant driving frequency. Sec- 
tions III and IV deal with theory of chirped Faraday 
waves, in inviscid (Sec. Ill) and low-viscosity (Sec. IV) 
fluids. Section V presents a brief discussion of our results. 
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II. WEAKLY NONLINEAR FARADAY WAVES 

A. Inviscid fluid 

To set the stage for a theory of chirped Faraday waves, 
we need to briefly review the theory of weakly nonlinear 
Faraday waves with a constant driving frequency. Con- 
sider a quasi-two-dimensional rectangular tank with a 
fluid of length I, width w and depth h, so that I ^> w. We 
assume that the elevation of the fluid, caused by the wave, 
depends only on the longitudinal coordinate x and time 
t, so that we have a quasi-two-dimensional flow in the 
xz plane (z is the vertical coordinate). The unperturbed 
level of the fluid is at z = 0. The vertical displacement 
of the vibrating tank is described by the equation 



(,{t) = ao cos(2u; t) . 



(1) 



We assume weak forcing, that is the vibration acceler- 
ation is much less than the gravity acceleration <?, and 
introduce a small dimcnsionless parameter e: 



^ 2 ao , 
e = < 1 . 



(2) 



In the limit of inviscid fluid the flow remains potential 
once it is potential at t = 0, and the external forces 
are potential [13, as is the case here. The assump- 
tion of a potential flow is also approximately valid in 
a low- viscosity fluid . We also assume that the wave- 
length of the standing wave is much larger than the cap- 
illary length of the fluid and neglect the capillary effects 
throughout the paper. The linear dispersion relation for 
the wave is w„ = gk n tanh(fc„/i), where to n is the natural 
frequency of the n-th mode, k n = 2ir/\ n = nir/l is the 
wave number of the n-th mode, and n = 1,2,.... 

The governing equations for the velocit y p otential 



v>t + ^(vl + v>t) + (g + Ov 



z = i) 

u 

=— ft, 



MM 




0, 


(3) 


0, 


(4) 


0, 


(5) 


0. 


(6) 



where indices denote partial derivatives. The Laplace's 
equation describes a potential flow of an incompress- 
ible fluid. Equations (@J and 10 are the Navier-Stokes 
equation and the kinematic boundary condition, respec- 
tively, evaluated at the free surface. Finally, Eq. © is 
the boundary condition for the vertical velocity compo- 
nent at the bottom of the tank. 

Let the vibration frequency be close to twice the nat- 
ural frequency of the primary mode n — 1, i.e. to ~ tui, 
so that this mode is excited via parametric resonance. In 
a weakly nonlinear regime it suffices to account for the 
excitation of only one higher order mode: the secon dary 
mode n — 2, which is enslaved to the primary mode |12| . 
Therefore, one should look for ip(x,z,t) and rj(x,t) in 



Eqs. OJl-® in the following form [TlIT^: 

cosh [ki{z + h)\ 



ip{x,z,t) = ip (t) + ipi(t)ipi(x) 



cosh(kih) 



r)(x,t) 



. . , , . cosh \k 2 (z+h)] 

MQMx) Zl m + 

cosh(fc 2 /i) 

T) l {t)<<]) l [x)+<q 2 {t)l]) 2 {x) + ... 



(J) 



where the eigenf unctions ip n (x) — V2 cos (k n x). The 
higher order terms will be neglected in the following. In 
the deep-water limit k n h^> 1, the linear dispersion rela- 
tion for the wave becomes ~ gk n . For this approx- 
imation to hold with an error less than 0.5%, it suffices 
to demand that h > I. 

The perturbation theory we are using employs the 
smallness of e. As will be seen later, this smallness 
implies a smallness of the wave amplitude compared 
with the wavelength, so that the dimensionless param- 
eter « = k\r\ is small. Expanding (p(x, z = rj,t) in the 
vicinity of the unperturbed surface z = in a power 
series in k, and substituting it and the second of Eqs. 
O into Eqs. (0} and JSJ, we obtain in the leading and 
sub-leading orders of k: 



m(t) 



V2 



m 



(8) 



(9) 



As we see, the next-order corrections (po and 772 are en- 
slaved to the primary mode, and their magnitudes are 
0(k?7i). In addition, we obtain a nonlinear differential 
equation of the second-order for the time-dependent am- 
plitude of the primary mode 771 (i) : 

m + \ k i (5vhi ~ Sufril) + ojj[l + 4e cos(2w t)] 771 = , 

(10) 

where we have used Eq. Q and kept terms up to C(k 3 ). 
Equation 110|> is a generalization of the linear Mathicu 
equation |20| . 

Now we employ the method of averaging .20] ■ We make 
an Ansatz 771 (t) = A\(t) cos [ivit + </>i(t)] and 771 (t) = 
— LUiAi(t) sin [uj\t + 4>i(t)] m Eq. (fTU|l and treat the am- 
plitude A\(t) and phase 4>i(t) as slow functions of time. 
(Being interested in the first-order equations with respect 
to g, one can omit higher temporal harmonics in 771 (£) 
|20ll21| .') Let 5 = lui—lu be the (small) detuning from the 
exact linear resonance. Then, for e <C 1 and \8\ <C Wi, 
we can perform averaging over the fast time ~ w-j -1 |20j . 
Introducing a new phase variable <f) — 5-t+<f>i, we obtain: 



Ai = euJiAi sin(2</>) , 
4> = £LO\ cos(20) 



k\A\uJi 



8. 



(11) 



The second term in the right side of the equation for <p 
describes the nonlinear frequency shift of the standing 
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wave. One can see that, as the wave amplitude grows, 
its frequency goes down j22|. This fact is important in 
the autoresonance excitation scheme introduced below. 
Rescaling time, amplitude and detuning, 



T — £U>it , B 



2^1 



-4, 



A 



8 

euii 



(12) 



we rewrite Eqs. <|11[) in a scaled form: 

B = £sin(2<£), 

4> = cos(20) - B 2 + . 



(13) 



where the time derivatives are taken with respect to the 
slow time r. When A < C(l), the typical value of B (for 
example the stable fixed point, see below) is 0(1). Going 
back to Eq. 1|12|) . we see that, in the dimensional units, 
the parameter kiAi ~ e 1 / 2 -C 1. As in the leading order 
k — kiAi, this validates our assumption that n <C 1. 

Equations H13|) describe weakly-nonlinear constant- 
frequency Faraday waves in the leading order in e. In 
the context of Faraday waves, Eqns. I|13(l were first ob- 
tained by Miles 0, though he derived them in a 
different way, working with the Lagrangian of the fluid. 
In the sub-lea ding order in e, additional nonlinear terms 
appear 

dill EH El, 

which will not be considered here. 
Equations (|13fl can be rewritten in a Hamiltonian 
form if we introduce the action and angle variables / = 
B 2 /2 and <p: 



I = -^ = 2/sin(20), 
H 

4> = — = cos(20) -21 + A. 
where the Hamilton's function is 



H(I, 



I cos(2(f>) - I 2 + A / . 



(14) 



(15) 



The fixed points of this dynamical system are deter- 
mined by the value of the scaled detuning A [23, |2^| : 

(a) A < — 1. No fixed points. 

(b) — 1 < A < 1. Three fixed points: an elliptic 
point [/» , 0*] = [(1 + A) /2, 0] and two saddle points 
[/»,</>*] = [0,±arccos(-A)/2]. 

(c) A > 1. Two fixed points: the same elliptic point 
[J*, 0*] = [(1 + A)/2, 0] as in case (b), and a saddle 
point [h,(j>*\ = [(A-1)/2,tt/2]- 

Figure n shows the phase plane (<j>,I) in the cases of < 
A < 1 and A > 1. The phase portrait is periodic in <j> 
with period tt. In the case of < A < 1, the separatrix is 
formed by the curve / = A + cos(2</>) and the straight line 
1 = 0. In the case of A > 1, the separatrix is formed by 
the curves / = [cos(20) + A — \/a\/2 and / = [cos(20) + 
A + v / a]/2, where a = cos 2 (2</>) + 2Acos(20) + 2A - 1. 
Notice that the maximum possible amplitude of phase- 
locked oscillations is achieved at a nonzero detuning from 




FIG. 1: The phase portrait of the inviscid constant-frequency 
system (1141 with detuning A = 0.5 (a) and A = 2 (b). Phase 
locking occurs inside the regions limited by the separatrix 
(denoted by the thick line). The saddle points are = 
[0, ± arccos(— A/2)] in case (a) and [1/2, ±7r/2] in case (b). 
The phase portrait is periodic in <f> with period n. 



the exact linear resonance, like in many other instances 
of nonlinear resonance. 

As the Hamilton's function (|15|) is a constant of mo- 
tion, the system is integrable. In particular, one can find 
the "nonlinear period" : the period of motion along a 
closed trajectory in the phase plane. Denoting the con- 
stant Hamilton's function as Hq, we obtain 



T n l — 2 



|[A + cos(20)] 2 -4i? } 



1/2 



(16) 



where <p± — ± arccos(2V Hq — A)/2. For a zero detuning, 
and initial conditions very close to the fixed point — 
1/2 and 4>* = (so that H — 1/4), we obtain T ~ n. 
This corresponds to small harmonic oscillations around 
the elliptic fixed point. In the physical units the period of 
small oscillations is ~ 7r/(ewi), that is much longer 
than the wave period. 



B. Low- viscosity fluid 

Taking into account a weak damping of the wave 
amounts to adding a linear damping term 27771 to the left 
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side of Eq. IjlOfl. where 7 is defined in terms of the rate 
of loss of mechanical energy due to dissipation, 0. The 
incorporation of only a linear damping term requires that 
j/u>i < l,so that the damping is treated perturbatively. 
The specific damping mechanisms which contribute to 
the value of damping rate 7 are the bulk viscosity [Isj . 
dissipation in the vicinity of the fixed walls |24j , dissipa- 
tion at the free surface (especially if contaminated) [24| 
and contact line damping |2f| (see, e.g. Ref. [23 for 
a review). In practice, one can interpret the damping 
rate as a phenomenological term, and determine it from 
a comparison with experiment. 



(a) A < — 1. No fixed points. 




FIG. 2: The phase portrait of the constant-frequency system 
1171 1 with a zero detuning and scaled damping rate V — 0.2. 
The thick lines mark the separatrices, which pass through 
the saddle points [_B*,0*] = [0,±7r/4]. The phase portrait is 
periodic in (j> with period 7r. 

Including the linear damping term in the first of Eqs. 
(|13fl . we obtain: 



B = Bsm(2(j>) - TB, 
<j> = cos(20) - B 2 + A . 



(17) 



where T = 7/(6^1) > is a dimensionless damping rate. 
It follows from the first of Eqs. (|1 T|) . that non-trivial 
fixed points £?* ^ can exist only when Y < 1, that is 
for a small enough viscosity. In this case, any trajectory 
on the phase plane of the system (except for trajectories 
with a zero measure) converges to a stable focus, see Fig. 
121 This is in contrast to the inviscid case, where start- 
ing from initial conditions outside the separatrix leaves 
the trajectory phase- unlocked. The fixed points [B*,^*] 
of Eqs. (|T7|l are determined by the values of A and V. 
Since Eqs. I|17(l [and Eqs. 1)13(1 ] are invariant under the 
transformation B — > — B, one needs to consider only fixed 
points with positive amplitudes. Let a = \/l — T 2 and 
£ = vl — A 2 . Let us also denote the critical damping 
rate 
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1/2 



(b) -1 < A < 0. For T > £, [0, arccos(-A)/2] is 
a stable node, and [0, — arccos(— A)/2] is a saddle 
point. For < T < £, [0, ± arccos(-A)/2] are 
two saddle points, and [(a + A) 1 / 2 , arcsin(r)/2] is 
a stable fixed point. For < T < T cr it is a stable 
focus, while for r cr < T < £ it is a stable node. 

(c) < A < 1. For T > 1, [0,arccos(-A)/2] is 
a stable node, and [0, — arccos(— A)/2] is a sad- 
dle point. For f < T < 1, [0, arccos(-A)/2] 
is a stable node, [0, — arccos(— A)/2] and [(A — 
er) 1//2 , 7r/2 — arcsin(r)/2] are two saddle points, and 
[(<7 + A) 1 / 2 , arcsin(r)/2] is a stable fixed point. 
For £ < T < r cr it is a stable focus, while for 
T cr < T < 1 it is a stable node. For < T < 
£, [0,± arccos(— A)/2] are two saddle points, and 
[{a + A) 1 / 2 , arcsin(L)/2] is a stable focus. 

(d) A > 1. For < T < 1, + A) 1 / 2 , arcsin(r)/2] is 
a stable fixed point. For < T < T cr it is a stable 
focus, while for T cr < V < 1 it is a stable node, and 
[(A - er) 1 / 2 , 7r/2 - arcsin(r)/2] is a saddle point. 

Figure [3] shows two characteristic values of the scaled 
damping T as functions of A. The first one is T cr from 
Eq. i|18fl • The second one is the maximum value of T for 
which a nontrivial stable fixed point still exists. For — 1 < 
A < this maximum value is equal to VT— A 2 , while for 
A > it is equal to 1. To conclude the brief review of the 
constant-frequency theory, we notice that the dependence 
of on L exhibits a pitchfork bifurcation. Figure 0] 
shows the bifurcation diagram in case (c) for A = 0.5. 
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FIG. 3: The solid line shows the critical scaled damping rate 
r cr as a function of the scaled detuning A [Eq. USU I. For 
r > r cr (the over-damped case) a stable node is obtained, 
while for F < r cr a stable focus is obtained. The dashed line 
shows the maximum value of F for which a nontrivial stable 
fixed point still exists. F cr is always below this maximum 
value. 



which will appear shortly. There are four different cases: 



5 



1.5 



1 

GO 

0.5 



0.5 V3/2 1 

FIG. 4: A combined bifurcation diagram of the constant- 
frequency system 1171 . Shown is the fixed-point amplitude B* 
versus T for A = 0.5. The solid line shows the stable focus or 
stable node (depending on F) , the dashed line shows the sta- 
ble node and the saddle point [as described in case (c)], and 
the dotted lines show the unstable fixed points (saddle points 
and unstable focus). The bifurcation occurs at V = yl — A 2 , 
where the stable node [B*,<f)*] — [0, (1/2) arccos(— A)] be- 
comes a saddle point. In the region yfi — A 2 < F < 1 both 
the stable focus, and the stable node exist. 



III. CHIRPED FARADAY WAVES IN AN 
INVISCID FLUID 

A. Governing equations, phase portrait, criteria 
and numerical examples 

Now let the vibration frequency be time-dependent 
(chirped): u) = u>(t). In general, the dimensionless pa- 
rameter e = uj 2 ao/g will also become time-dependent. 
For simplicity, we shall assume that ao also varies in time 
so that e = const . Our objective is to keep a Faraday 
wave close to resonance in spite of its nonlinear frequency 
shift, so as to achieve a persistent growth of the wave am- 
plitude. Like in other autoresonance schemes, the exact 
form of the function u>(t) is unimportant if this function 
satisfies three criteria: 

1. The chirp sign coincides with the sign of the non- 
linear frequency shift of the wave. For the standing 
Faraday waves ui(t) should decrease for the wave 
amplitude to increase. 

2. The frequency chirp rate must be sufficiently small, 
so that the phase portrait of the system evolves 
adiabatically: \w{t)\T n i <C where T n i is the 
characteristic nonlinear period, see Eq. I|lu|) . 

3. The dynamic frequency mismatch, which we de- 
fine as the absolute value of the increment of the 
vibration frequency during one nonlinear period, 
\u)(t + T n i) — u)(t)\, should be small compared with 
the inverse nonlinear period, T~, . In physical 
units, this yields ^/(euii) 2 <C 1. 



Criteria 1 and 2 have appeared in previous works on 
autoresonance 0], while criterion 3 is new. We shall see 
shortly that, in the problem of parametric autoresonance, 
criterion 3 is more restrictive than criterion 2. 

The derivation of the equation of motion for the pri- 
mary mode amplitude, for a slowly time-dependent driv- 
ing frequency, goes along the same lines as in the case of 
a constant driving frequency. The resulting equation is 
(compare with Eq. (|10fl ]: 

Vi + 2 k K 5 rtvi - 3u> 2 r)l)+u> 2 [1 + 4ecos(2$(i))] m = , 

(19) 

where $(t) = f Q u>(t') dt' . If e is small, and the chirp 
rate is slow on the time scale of the wave period, one can 
again use, close to the parametric resonance, the method 
of averaging [?J . For concreteness, we assume in this work 
a constant chirp rate \i: 

uj{t) = u> 1 - fit , (20) 

so that Q(t) — uj\t — /j,t 2 /2. Introducing a scaled chirp 
rate m = u/(ujie) 2 , and the same scaled time r and am- 
plitude B as before [see Eq. (|12|) ]. we arrive at the fol- 
lowing scaled equations: 

B = B$m(2(t>), 

<j) = cos(2(f)) - B 2 + mr, (21) 

where now <f>(t) = fit 2 /2 + 4>i(t), and the differentiation 
is done with respect to r. In the action-angle variables 
we obtain: 

I = 2/sin(20), 

4> = cos(24>)-2I + mT. (22) 

Once B(t) and 4>(t) are found, one can immediately re- 
construct the standing way profile t)(x, t) by using the 
Ansatz r)i(t) = Ai(t) cos [u>it + 4>i{t)] and the "enslaving 
relation" JSJl in the second of equations J7J). Therefore, 
in the rest of the paper we shall focus on Eqs. (1221 which 
coincide, up to notation, with those obtained by Khain 
and Meerson 0, who investigated parametric autoreso- 
nance in a nonlinear oscillator. Equations similar to (1221) 
also appear in the problem of the second-harmonic au- 
toresonance in an externally driven oscillator |27| . 

The Hamilton's function of the system (1221) is time 
dependent: 

H{I, <f>) = I cos(20) -I 2 + mrl , (23) 

so H is not a constant of motion anymore. In the follow- 
ing we shall use t for the slow time r. 

Figure [S] shows an example of parametric autores- 
onance: a persistent phase locking and a systematic 
growth of I with time, with some oscillations on top of 
the systematic growth. Here the scaled chirp rate m is 
less than some critical value m cr for these initial condi- 
tions. Figure illustrates breakdown of autoresonance, 
observed when m > m cr . As in other instances of au- 
toresonance, a theory of parametric autoresonance ap- 
peals to the constant- frequency theory. Comparing Eqs. 
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FIG. 5: Parametric autoresonance in an inviscid fluid. Shown 
are the action 7 (a) and the phase cf> (b) as functions of time. 
The parameters are m = 0.2, 7(0) = 0.4, <f>(0) = 0, and a zero 
initial detuning. The solid lines shows numerical solutions. 
The dashed lines show the trends 7* (t) , given byEq. @U (a) 
and <j>*(t), given by Eq. H25H (b). 

(|22[) and (|14|) . one can see that the term mt plays the 
role of an effective (time-dependent) detuning. There- 
fore, when the chirp rate is small, m <ti 1, the phase 
portrait of the system almost coincides with that of the 
autonomous equations 114|) (see Fig. 1), except that now 
it changes with time according to the current value of the 
effective detuning. The change of the phase portrait is 
adiabatically slow, except at times t ~ —1/m and 1/m 
(corresponding to A = — 1 and 1, respectively), when bi- 
furcations occur. One consequence of the adiabatic evo- 
lution is that Eqs. I122[I have "quasi-fixed" points. The 
most important stable quasi-fixed point [I*(t),<j)*(t)} can 
be found by assuming that <p*(t) <C 1, and that it varies 
with time slowly. Then the second of Eqs. yields, in 
the leading order, 

I* - (1 + mt) . (24) 
Substituting this into the first of Eqs. Q22JI. we obtain 



The stable quasi- fixed point, or trends (|24l) and (|25H . 
previously found by Khain and Meerson Q (see also 



Ref. [23), are the essence of parametric autoresonance. 
Shown in Fig. [S] are I(t) and <f>(t) found numerically, 
and the trends and The trend J2U cor- 

responds to a steady growth of the wave amplitude: 
B*(t) = [21* (i)] 1 / 2 ~ (mt + I) 1 / 2 . The important phase 
trend Q25[l was overlooked in Ref. j^]. Notice that, at 
scaled time t > 1, the phase trend 0* ~ 1/(4^) becomes 
independent of the chirp rate m. Importantly, for the 
expressions (I24|l and (|25|l to be valid, one can either de- 
mand m <C 1, or go to long times: t 3> 1. Therefore, the 
stable quasi-fixed point keeps its meaning, at long times, 
even at finite (non-small) m. 

We found that, surprisingly, unstable quasi- fixed points 
of the chirped system also play an important role in 
the dynamics. The unstable points are analogs of the 
constant-frequency saddle points discussed in the pre- 
vious section [see the text following Eq. I|15[l]. To 
find the locations of the unstable quasi-fixed points 
in the leading order, one can simply replace the de- 
tuning A by mt. Therefore, on the time interval 
< t < 1/m, there are two saddle quasi-fixed points 
[/*,</>*] — [0, ±(1/2) arccos(— mt)]. These points dis- 
appear at t ~ 1/m, and a new saddle points appears: 
[/*,</>*] — [(mt — l)/2,7r/2)]. These expressions (includ- 
ing the boundaries of the corresponding time intervals) 
are valid in the leading order in m < 1. Higher-order 
corrections can be also calculated. 

Now we are in a position to discuss criteria 2 and 3 
for parametric autoresonance in this system. For a con- 
stant chirp rate fj, [see Eq. (|2U|l ]. criterion 2 can be writ- 
ten, in the physical units, as [i(euj)~ l < u, or // < eu> 2 . 
Now, the dynamic frequency mismatch, acquired by the 
chirped system during time T n i, can be estimated as 
fiT n i ~ [i/(eu)). Criterion 3 demands that this quantity 
be small compared to ~ em, which yields fx < e 2 uj 2 . 
As e is small, criterion 3 is more restrictive than crite- 
rion 2. In the scaled variables, criterion 3 has the form of 
m < 1, as can be expected from the form of scaled Eqs. 
(|22|l . The inequalities here are written up to numerical 




FIG. 6: Breakdown of parametric autoresonance in an inviscid 
fluid in the case of m = 5.5 > m cr . The parameters are m = 
5.5, 7(0) = 0.6, 0(0) = 0, and a zero initial detuning. Shown 
are the action I(t) and the phase cj>(t) found numerically. 
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factors which depend on the initial conditions, see below. 

One more convenient description of the chirped system 
can be achieved if we rewrite Eqs. (|22|l as a second order 
equation for the phase: 

+ 2mt sin(20) + sin(40) - m = , (26) 

or 

where we have introduced a time-dependent potential 

V(<p,t) = --cos(4</>) - mi cos(20) - m0. (28) 

This suggests new canonical variables and u = 
cos(2cf>) — 21 + mt, so that in the new time-dependent 
Hamiltonian, H ((/>, u, t) = u 2 /2 + V(<fi, t), there is a clear 
separation between the potential energy and the kinetic 
energy. The new Hamiltonian describes a "particle" of 
a unit mass and velocity u = <f>, moving in a time- 
dependent potential V. This picture is useful for a quali- 
tative analysis of the dynamics of the "particle" when to 
is small, so the potential slowly varies in time, see Fig. [7| 
In the variables u, <j) the stable quasi-fixed point becomes 
(approximately) [0, (m/4)(l + mt)~ 1 ], while the saddle 
points are [0, ±(1/2) arccos(— mt)] at < t < 1/m, and 
[0, ~ ± 7r/2] at t > 1/m. We shall see shortly that each of 
the unstable points [0, (1/2) arccos(— mt)] and [0, ~ 7r/2] 
plays an important role in this system. 




FIG. 7: The time-dependent potential V(<p,t), Eq. is 
shown as a function of 4> f° r several consecutive times start- 
ing from t = 0. The scaled chirped rate m = 0.1. The po- 
tential well at <f> ~ 4>*{t) deepens with time starting from 
t = 0. At < t < 1/m there is a potential barrier at 
4> ~ (1/2) arccos(— mi) (denoted by the circle), which dis- 
appears at t ~ 1/m. At t ~ 1/m a new potential barrier 
appears close to cj> = tt/2 and heightens with time (denoted 
by the asterisk). 

Let us consider two typical cases of parametric au- 
toresonant excitation of a Faraday wave. In the first 
case one first excites the wave at a constant frequency, 



so that the initial values of the action and phase are 
in the vicinity of the stable fixed point. Then, upon 
slowly reducing the driving frequency, one keeps the 
phase locked, as our "particle" oscillates in a potential 
well which slowly deepens with time, see Fig. [7] In 
the second case one starts the autoresonant driving from 
an almost zero wave amplitude. Here the saddle point 
[/*, <f>*] ~ [0, (1/2) arccos(— mt)] plays an important role. 
Indeed, the stable manifold of this quasi-fixed point is 
along the <f> axis. Therefore, no matter what the initial 
phase is, the phase approaches, on a time scale 0(1), 
the saddle point. The unstable manifold of this saddle 
point is along the / axis, so I(t) will grow with time. 
Still, if I(t = 0) is small enough, I(t) remains small 
during this time scale 0(1). Therefore, phase locking 
is always achieved at this stage, so the time interval 
< t < 1/m can be called the "trapping stage". Later 
on I(t) grows significantly but, as we found numerically, 
the "particles" remain inside the (slowly expanding) sep- 
aratrix / = cos(20) + mt. As a result, the phase starts 
to perform large-amplitude oscillation around the stable 
quasi-fixed point, and phase locking persists. 

One more alternative description of the system of equa- 
tions is in terms of the complex amplitude i/)(t) = 
B(t)exp[i<f)(t)]: 

iipt + i>* - (H 2 - mt)ip = 0, (29) 

where the subscript t denotes differentiation with respect 
to the slow time. The long-time behavior of can be 
obtained by looking at the asymptotic solutions of Eq. 
(|29|l at t — > oo. For a solution such that \ip\ grows with 
time like a power law, the leading terms are those in the 
parentheses. This immediately yields 

l^l-fmi) 1 / 2 , (30) 

which corresponds to the leading term (when mt 1) 
of Eq. H24f> . and describes a phase- locked wave [here <j> 
stays close to zero, see Eq. (|25|) ]. On the contrary, if 
\tp(t)\ remains bounded and small, the first term of Eq. 
(I29|l is balanced by the last one, and we obtain 

iP(t) = Voexp (imt 2 /2) , (31) 

where ipo = \ipo\ e x p(^o) = const. This solution corre- 
sponds to an unlocked phase (f>(t) = (f>o + mt 2 /2 and a 
constant amplitude \ipo\- Of course, the phase of the wave 
<j>i(t) = <fi(t) — mt 2 /2, which is defined by the Ansatz 
r)i(t) = Ai(t) cos [u}\t + 4>i(t)] (where t is the physical 
time), stays constant in this regime, and is equal to <po. 

We determined numerically, for several typical classes 
of initial conditions, the critical value of to, m = m cr , 
which separates the phase locking regime from the phase 
unlocking regime [28|]. At a fixed initial phase <f>(t — 0) = 
0, m cr grows with the initial amplitude B(t — 0), at least 
until the scaled amplitude becomes of order unity, see 
Fig. El This result agrees with those obtained by Fajans 
et. al. (see Fig. 3 in Ref. |27|) who presented them 
in terms of e cr versus u. When starting from the stable 
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FIG. 8: The critical chirp rate a function of the initial 

amplitude B(0) for (j>(0) = 0. 



quasi-fixed point: I(t = 0) = 1/2 and 4>{t = 0) = 0, we 
found that m cr ~ 4.963. 

We also found TTl/Qf clS 0i function of the initial phase, 
for a given, and very small, initial amplitude, see Fig. [5] 
This dependence is relatively weak. The largest m cr is 
obtained for tj> = 7r/4, the smallest one for <f> = —ir/4. 




FIG. 10: Autoresonance and its breakdown for m very close 
to m cr ~ 4.963. The initial conditions are (f>(t = 0) = 
and I(t — 0) — 1/2. The dashed lines show the phase-locked 
solutions for the phase and action for m just below m cr , while 
the solid lines show the phase-unlocked solutions for m just 
above m cr - 



unlocking). Instead, the "particle" lands, at t = oo, on 
the peak of the potential at <f> — tt/2. In the following we 
shall find the asymptotic form of this special trajectory. 




-7172 



-71/4 , 



FIG. 9: The critical chirp rate a function of the initial 

phase 4>(0), for a very small initial amplitude -B(O) = 10 -6 . 

What is the signature of the special case to — m cr ? At 
to < m cr the phase oscillates in time. As to approaches 
m cr from below, the onset of the phase oscillations is de- 
layed more and more, see Fig. 1101 Now, at m > m cr the 
phase initially changes slowly and then rapidly escapes to 
infinity. As m approaches m cr from above, the point of 
rapid departure of the phase is delayed more and more, 
as shown in Fig. ^] This suggests that, in the spe- 
cial case to — m cr the phase neither oscillates around a 
trend, nor departs from it. Instead, the phase monotoni- 
cally approaches tt/2. Looking at the effective potential, 
shown in Fig. Q one realizes that, in the special case 
m = m cri our "particle" neither oscillates in the poten- 
tial well (which would correspond to phase locking), nor 
escapes from the well (which would correspond to phase 



B. Perturbative solutions 

In this subsection we present three perturbative ana- 
lytic solutions which illustrate the basic features of para- 
metric autoresonance: persistent resonant growth, cap- 
ture into resonance and the limiting trajectory which sep- 
arates between phase locking and unlocking. In each of 
the three cases, a local analysis around one of the quasi- 
fixed points of the system is required. 



1. In the vicinity of the stable quasi-fixed point 

Let us linearize Eq. (|26|l in the vicinity of (/) = 0. As 
discussed above, this requires either a small chirp rate, 
to <C 1, or a long time, 4 > 1. We obtain: 



<fi + 4(1 + mt)(f> = to . 



We look for the solution in the form 

TO 



4(1 + rat) 



84>(t) 



(32) 



(33) 



where the first term is the phase trend (|25|l . Neglecting 
higher-order terms, we arrive at the Airy equation 



5(f) + 4(1 + mt)5(j} = , 
whose general solution is [29] ]: 



(34) 



C 2 Bi 



2/3 



(1 + mi) 



2/3 



(1 + mi) 



(35) 
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Here Ai(r) and Bi(r) are the Airy functions of the first 
and second kind, and C\ and Ci are constants depend- 
ing on the initial conditions. Using the large-argument 
expansion of Ai(—z) and Bi(—z) [29j, we obtain: 



A 



4(1 + mi) (l + mt)V4 



4(l+mi) 3 / 2 



3m 



(36) 

where A and £ are new constants. The action I(t) can 
be found using the second of Eqs. 



J(t) 



mi + 1 



A(l + mi) 1/4 < 



4(1 +mf) 3 / 2 



3m 



(37) 

The first terms in Eqs. (|36|l and l|37|) are the systematic 
trends (|2EJ> and ifiHjl. The solutions (jSEJ and coin- 
cide, up to notation, with the WKB-solutions obtained 
in Ref. Q. Fig. ^l s h ows excellent agreement between 
the analytical solutions (|3l)|) and l|3TI) and numerical so- 
lutions. 



2. Capture into resonance 

Now we consider the case of driving a Faraday wave 
starting from a very small amplitude and a large ini- 



0.05 



-0.05 




FIG. 11: The phase <j>(t) (a) and the action I(t) (b) as func- 
tions of time for m = 0.2, 0(0) = and 7(0) = 1/2. The 
asterisks mark the asymptotic solutions given by Eqs. (1361 1 
(a) and 1371 (b). The solid lines mark the numerical solutions. 



tial detuning, that is far from resonance. Here one is 
interested in the capture into resonance. In the case of 
external autoresonance this phenomenon has been exten- 
sively studied by Friedland U • In the case of parametric 
autoresonance this phenomenon has not been addressed, 
although equations similar to (|22[1 were analyzed in Ref. 
|27| in the context of the second-harmonic autoresonance 
in an externally driven oscillator. 

We start, in Eqs. l|2Tjl. at a large negative time to < 0, 
and assume that the initial detuning mtQ is very large 
in the absolute value, while the initial scaled amplitude 
B(t = to) is much less than unity. As long as B(t) -C 1, 
one can neglect the £? 2 -term in the second of Eqs. 121|) : 



cb — cos(20) + mi . 



(38) 



This equation describes two distinct stages of the dy- 
namics. In the first stage i is large and negative, and 
\mt\ ^> 1. Therefore, the term mt is dominant, so 27] 



0~0(i o ) + -(i-i o ) 2 



(39) 



During this pre-locking stage, the phase varies rapidly. 
The second stage occurs roughly at — 1/m < t < 1/m. 
Here the two terms on the right hand side of Eq. i|38|) 
are comparable. As m <C 1, the duration of this stage is 
long compared with unity, and the phase approaches, on 
a time scale 0(1), the unstable quasi- fixed point 



arccos(— mi) 



4[1 - (mt) 2 } 1 / 2 



(40) 



where we have included the next-order correction in m. 
As noted previously, this quasi-fixed saddle point ceases 
to exist at t > 1/m; the correction term in Eq. (|40[1 is 
invalid too close to t = 1/m. We call this regime "linear 
phase locking" stage. 

Let us obtain the dependence of B on time in the 
pre-locking and linear phase locking stages. In the pre- 
locking stage 4>(t) is given by Eq. I|39|) . Then the first of 
Eqs. yields 

B(t) ~ B(to) exp |y sin [m (s - t ) 2 + 20(t o )] ds 

(41) 

(the integral can be expressed through the Fresnel inte- 
gral). At this stage B(t) oscillates rapidly, because of the 
rapid and monotonic change of the phase. At the linear 
phase locking stage, <p(t) ~ (j)*(t) as given by Eq. JJOJ. 
Then, neglecting higher-order terms in m, we arrive at 
the equation 



which yields 



B ~ £?2 exp 



B ~ By/1 - (mt) 2 

arcsin(mt) 



ty/i - (mt) 2 + 



(42) 



(43) 



where B2 is a constant determined by the initial condi- 
tions. 
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To find the trend, we neglect the term Sep. Therefore, at 
t ^> 1, we obtain Sep ~ 1/ (4i), so 



7T 1 

2 ~ it 



(45) 



Using the second of Eqs. I|22() . we obtain the respective 
trend of I(t): 



. . mt 1 



(46) 



At this stage we notice that Eqs. (|45|l and (|46() describe, 
at t 3> 1, one of the unstable (saddle) quasi- fixed points 
of the system: the one that appears close to t = 1/m. 
This again shows that adiabaticity holds not only at m <C 
1, but also at t 3> 1. 



71/2 




FIG. 12: The dynamics of the system in the linear phase 
locking stage —1/m < t < 1/m. The phase <p(t) (a) and 
amplitude B(t) (b), found numerically, are shown by solid 
lines. The quasi- fixed saddle point (j>*(t) [Eq. I40H 1 (a) and 
the amplitude [Eq. (0HI] (b) are denoted by dashed lines. The 
initial conditions are to = —400, 4>(~to) = 0, and B(to) = 10 -6 , 
the scaled chirp rate m = 0.01. It can be seen that, as B 
becomes comparable to 1, the linear theory breaks down. 

Equation l|43|) breaks down when the earliest of the 
two events occurs: t becomes larger than 1/m, or B(t) 
becomes comparable to unity, so that one cannot neglect 
the £? 2 -term in the second of Eqs. (|21|l . The phase 
starts to oscillate, while the amplitude both oscillates 
and grows, see Fig. ^1 and the system enters the au- 
toresonance regime. 

Similar results are observed when starting from a small 
amplitude and a zero initial detuning (that is, in exact 
linear resonance). Therefore, when the initial amplitude 
is very small, and m« 1, phase locking is very robust. 

3. Critical chirp rate and the limiting trajectory 

As we have seen numerically, when m approaches the 
critical value m cr , the phase <fi approaches 7r/2 at t > 1. 
Let us assume that, at a time to, (f> is already in the 
vicinity oiir/2: (j) = tt/2 — 8<j), where 84> <C 1. Linearizing 
Eq. (123), we obtain 

S"4> - 4:5<j)(mt - 1) = -m . (44) 



2 t 4 

FIG. 13: Numerical solutions for the phase 4>(t), for three 
different initial conditions, versus the universal asymptotics 
145H . denoted by the dotted line. The solid, dashed-dotted 
and dashed lines show the phase dynamics when starting from 
</>(0) = (here m ~ m cr — 4.963), 1/4 (here m ~ m cr ~ 
6.237), and 1/2 (here m ~ m cr ~ 7.094), respectively. 1(0) = 
1/2 in all three cases. 




2 4,6 8 10 



FIG. 14: Shown is I(t) for m just below m cr ~ 4.963. The 
solid line shows the numerical solution, starting from 1(0) = 
1/2 and (^>(0) = 0. The dashed line shows the trend l|46|l . 

Now consider small deviations from the unstable trends 
g5J and gnj>- Putting <j) = tt/2 -l/(4t) -£</>(*), we again 
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arrive at the Airy equation for 5(j)(t). Its general solution 
is 8<P{t) = dAi^mf/H] + C 2 Bi[(Am) 1 / 3 t}, where C x 
and C2 are determined by the initial conditions. Assum- 
ing m}/ 3 t 3> 1, we employ the asymptotics of the Airy 
functions |29|: 



Ai(z) 



exp [-(2/3)z 3 / 2 ] 



Bi{z) 



exp [(2/3) 



,3/21 



,1/4 



(47) 

at z > 1. The exponentially decaying term is negligi- 
ble. The exponentially growing term describes instability 
of the special trajectory with respect to small perturba- 
tions. This instability occurs both at to < m cr (when 
the phase leaves the vicinities of the saddle point and 
goes to the vicinity of the stable quasi-fixed point), and 
at to > m cr (when the phase locking terminates). For 
the special trajectory, obtained for to = to c ,., the coeffi- 
cient C2 must vanish, which brings us back to Eqs. H45|) 
and 1461) . where we must put m — m cr . Now it is clear 
that the asymptotic behavior of the 4>* (t) at to = m cr is 
independent of m cr and, therefore, on the initial condi- 
tions. On the contrary, the asymptotic behavior of /*(£) 
at to = m cr does depend on m cr and, therefore, on the 
initial conditions. Unfortunately, the local analysis does 
not enable one to find the value of m cr analytically. 

Figures ED an d 101 show the behavior of the system 
when to is very close to (just below) m cr . Figurc lTSl shows 
that, after a time 0(1), the phase follows the trend (|45ll . 
independently of the initial conditions (and of the value 
of m cr ). FigurelTHshows that, after a time of 0(1), good 
agreement between I(t), found numerically, and the trend 
(|46|l holds until the time when the "particle" leaves the 
vicinity of the unstable point and transfers to the vicinity 
of the stable point. 



IV. CHIRPED FARADAY WAVES IN A 
LOW- VISCOSITY FLUID 

We now account for a small viscosity and briefly de- 
scribe the dynamics of the phase and amplitude of the 
primary mode in the case of a slowly time-dependent vi- 
bration frequency. The weakly nonlinear governing equa- 
tions are obtained by analogy with Eqs. I|17fl : 



B = 5sin(2eV) -TB, 
4> = cos(20) - B? +mt, 



(48) 



where t is the slow time as before. For T < 1, there exists 
a stable quasi-fixed point which describes autoresonant 
excitation of the wave: 



(i-r 2 ) 1/2 + mt 



1/2 



1 . TO 

— arcsm(r) H ; - ; „ % ; 

2 K J 4[i-r 2 + TOi(i-r 2 ) 1 /2] 



When r < 1, Eqs. (j2<3) become 



B* 



Vl + mt , 0* ~ - + — — -r 

2 4(1 + mt) 



(49) 



(50) 




FIG. 15: Parametric autoresonance in the presence of viscos- 
ity. Shown is a projection of a three-dimensional trajectory in 
the space of <f>, B and t, upon the {(j>, £?)-plane. The solid line 
shows a numerical solution, the dashed line shows the trends 
148H . After several nonlinear cycles, the trajectory converges 
into a (slowly time-dependent) stable focus, which moves up- 
wards. The parameters are F = 0.2, m = 0.1, B(0) — 0.8, 
<f>(0) = 0.1, and a zero initial detuning. 



Figure ^| shows a projection on the (<fi, B) plane of a 
three-dimensional trajectory in the space of 0, B and t. 
One can see phase locking and a steady growth of the 
wave amplitude with time. 

Figures II til and 1171 show two different regimes of au- 
toresonant growth in the dissipative system. Shown in 
the two figures are the primary mode amplitude ver- 
sus time for the scaled damping rates Y < T cr and 
r > r cr , respectively. As the initial detuning is zero, 
T cr = 2/V5 = 0.89... here, see Eq. Figure 
[TBI shows decaying oscillations on top of the amplitude 
growth, given by the first of Eqs. (|49|l . Figure IT71 shows 
a non-oscillatory regime of the amplitude growth. Fig- 
ure !18l shows breakdown of autoresonance when the chirp 
rate m exceeds a critical value. 
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FIG. 16: Parametric autoresonance in the presence of viscos- 
ity. Shown is the amplitude B{i) versus time for a subcritical 
damping. The solid line denotes a numerical solution, while 
the dashed line shows the trend B,(t), see Eq. 14911 . The 
parameters are the same as in Fig. 1151 
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FIG. 17: Parametric autoresonance in the presence of viscos- 
ity. Shown is the amplitude B(t) versus time for a supercriti- 
cal damping. The solid line denotes the amplitude B(t) found 
numerically, and the dashed line denotes its trend B*(t), see 
Eq. Pty. The parameters are T = 0.9, m = 0.1, B(0) = 1, 
4>(0) — 0.25, and a zero initial detuning. 
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FIG. 18: Breakdown of autoresonance in the presence of vis- 
cosity. Shown are the amplitude B(t) and the phase </>(i), for 
m ~ 4.4 (just above m cr ) for the initial conditions B(0) = 1 
and 0(0) = 0.05, and for T — 0.1. One can see that, af- 
ter a strong transient excitation, the phase unlocks, and the 
amplitude decays. 



Similarly to the inviscid case, one can rewrite the gov- 
erning equations l|4*5|) as a single equation for the complex 
amplitude, tp = B exp(z0): 



iipt + V - (M 2 -mt- iT)ip = . 



(51) 



Assuming a solution growing in time (that is, phase 
locked solution), we obtain, at t 3> 1, B = \tl>\ ~ (mt) 1 / 2 , 
as before. 

Equations l|48|) can be also rewritten as a single second- 
order equation for (f>(t): 

4> + 2T<p + sin(40) + 2mi sin(20) 

- 2r cos(2</>) = 2Fmi + m (52) 

[compare it with Eq. (|26|) ]. This equation is convenient 
for a perturbative treatment in the vicinity of the stable 




FIG. 19: The phase <j>(t) (a) and the amplitude B(t) (b) as 
functions of time in the presence of small viscosity. The pa- 
rameters are m = 0.1, T — 0.1. The initial conditions are 
0(0) = 0.05 and B(0) = 0.95. The asterisks mark the analyt- 
ical solutions, the solid lines show numerical solutions. 



quasi- fixed point. For T -C 1, we can linearize Eq. I] ">-?!) 
around d> = 0: 



$ + 2f> + 4(1 + mt)(j> = 211(1 + mt) + m . 



(53) 



Now we substitute <f)(t) = (j)*(t) + 5<fi(t), where 0* is given 
by the second of Eqs. I|5U|) . and obtain a linear equation 



2TS(j) + 4(1 + mt)4> = 0. 



(54) 



Its approximate solution, at m <C 1 and T <C 1, can be 
written as 



Ae~ rt 
(l + mt) 1 / i 



sin 



4(1 + mt) 3 / 2 



3m 



(55) 



where A and £ are constants depending on the initial 
conditions. The respective solution for B(t) is 



B ~ (1 + mt) 



1/2 



Ae~ rt 
(1 + mi)V4 



4(1 + mtf/ 2 



3m 



cos h £ 

. ( ^ 6) 

These solutions are simple extensions of the non- viscous 
perturbative solutions. A comparison with numerical so- 
lutions is shown in Fig. 1191 and excellent agreement is 
observed. 
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In the viscous case the critical value m cr of the scaled 
chirp rate m = /z/(exo>i) 2 depends on the scaled damping 
rate of the wave T = ^/(euj-l). Figure I2T11 shows this de- 
pendence, which we found numerically when starting at 
t = from the stable fixed point of Eqs. lfPT|) with A = 0, 
that is from £?* = (1 — T 2 ) 1 / 4 and 0* = (1/2) arcsin(r). 
One can see that, as V increases, the critical chirp rate 
goes down monotonically. As T approaches 1, m cr goes 
to 0. Notice that, once we return to the physical (dimen- 
sional) critical chirp rate /i cr and the wave damping rate 
7, the dependence of [i cr on e, at fixed 7, is not a power 
law. 




r=Y/eco 

FIG. 20: The critical scaled chirp rate m cr = /i cr /(£ 2 ijjf) 
versus the scaled damping rate of the wave F = y/(eui), as 
described by Eqs. I|48|l . The initial amplitude and phase 
correspond to the stable fixed point of Eqs. (1171 with A = 0. 
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FIG. 21: B(t) (the thick solid line) and 771 (the thin solid 
line) versus time. The parameters are F = 0.006, m — 0.18, 
e = 0.006, 771(0) = B(0) = 0.9 and 7^(0) = <t>(0) = 0. 

Finally, we tested the accuracy of our reduced equa- 
tions |@SJ|. We compared numerical solutions of Eqs. 
(|48|l with numerical solutions of the unreduced equation 
of motion 119(1 for the primary mode amplitude ?7i(i), 
with a linear damping term added. Rescaling the time 
t = u>it, and the amplitude fji = (fci7yi)/(2y / e), one can 



rewrite the unreduced equation as 

•2 n 

fj 1 + 2eTfi 1 + lOsrj^i — 6er7j + 

+ fn [1 + 4e cos (2t - meV)] = . (57) 

A typical example of this comparison is shown in Fig. 
1211 and a fairly good agreement between the envelope of 
r)i(i) and the amplitude B{t) is observed. 

V. DISCUSSION 

This paper presents a theory of weakly nonlinear 
standing gravity waves, parametrically excited by weak 
vertical vibrations with a down-chirped vibration fre- 
quency. We have shown that autoresonancc phase locking 
and a steadily growing wave amplitude can be achieved 
despite the nonlinear frequency shift of the wave. For 
typical initial conditions we have found the critical chirp 
rate, above which autoresonance breaks down. When 
starting from a very small wave amplitude, and slowly 
passing through resonance, phase locking always occurs. 
We have obtained approximate analytical expressions for 
the time-dependent wave profile in different regimes. We 
have demonstrated that each of the three quasi-fixed 
points of the reduced dynamic equation, describing the 
primary mode, plays an important role in the dynamics 
of the system and/or in determining the critical chirp 
rate. 

Parametric autoresonance in Faraday waves is a robust 
phenomenon, and its experimental observation should 
not be difficult. To test our theory, one should use a 
quasi-two-dimensional tank with a low-viscosity liquid, 
and perform measurements of the standing wave eleva- 
tion as a function of time. The long-term wave amplitude 
trend [see the first of Eqs. (|49|) ]. and the critical chirp 
rate (see Fig. 20), give examples of quantitative predic- 
tions of the theory that can be tested in experiment. One 
such experiment is presently under way [30j. 

The applicability of the weakly nonlinear theory pre- 
sented in this work is limited to weak forcing and weak 
damping. Our analysis neglected hi ghe r order terms in 
e, which include a nonlinear forcing [l3|. a cubic damp- 
ing 0, 0| , a small correction to the linear detun- 
ing/frequency shift, a quintic conservative term [l^, etc. 
These terms can be included in the theory of autores- 
onance in order to achieve a better accuracy. Impor- 
tantly, once these higher-order terms are added to the 
constant-frequency model [Eq. (|17|) ]. the non-trivial sta- 
ble fixed point will exist onl y up to a certain value of the 
frequency detuning [lj, [lj, Ha, E3 ■ Therefore, in exper- 
iment, the autoresonant growth of the wave is expected 
to terminate when the (time-dependent) frequency de- 
tuning comes close to the maximum value, for which the 
non-trivial stable fixed point in the underlying constant- 
frequency model still exists. The breakdown of autoreso- 
nance is expected to occur at an amplitude smaller than 
the wave-breaking amplitude, which is comparable to the 
wavelength |3l) . 
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